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The problem of electron-proton scattering is handed over both the elastic and inelastic scattering. 
Two models are presented in this sense. The first, depends on the multi photon exchange ladder 
diagram, where the transition matrix is expanded in multi steps form. The second model uses 
the multi peripheral mechanism developed for the electromagnetic field. It allows the particle 
production in the inelastic scattering processes. An iterative procedure is found and inserted in 
a Monte Carlo program to reproduce the differential cross section of the reaction. The comparison 
with the experimental data shows bid fair in most cases. 



- - I. INTRODUCTION 

o ■ 
o ■ 

O ! Recent experiments for electron-proton (ep) scattering, SLAC-E-136 [1], DESY-HERA-HE [2], SLAC-E-133 [3] 
04 ■ and others [4], played important role in probing the nucleon structure and to reveal the dynamic mechanism of the 
^ [ electron interactions inside the nucleon bag. Many trials have been done before by a wide variety of empirical and 

theoretical models, ranging from form factor scaling [5], and vector meson dominance [6] to quark- parton models [7] 
'"J and perturbative QCD [8]. An eikonal optical picture [9-10] was used, based on multiple scattering of the incident 
CO , electron with the constituent valance quarks of the target hadron, assuming different forms of the (electron-quark) 

binary wave functions. 

[ In this article we shall deal with the problem of elastic and inelastic collisions of the electron with proton from 
a different point of view. In the following, we present two pictures for the (ep) scattering. The first is a scattering 
^ ^ in view of multi-photon exchange mechanism (MPEM) which is very convenient for the elastic scattering in a wide 
. range of momentum transfer square Q^. On the other hand, the multi-peripheral model (MPM) is relevant to the 
^-H inelastic and deep inelastic scattering and so works properly for the problems of multi particle production. The article 
, structure goes as follows. In section 2, we present the postulates of the MPEM, the formulation and the scenario 
■ of the MPM are given in section 3. The Monte Carlo generators are summarized in section 4. Finally results and 
] discussion are given in section 5, followed by conclusive remarks. 



II. THE MULTI PHOTON EXCHANGE MODEL (MPEM) 



D , 

' It is assumed that elastic scattering of electrons on proton proceeds via multi step processes represented by ladder 

^ I diagrams [11]. We proceed on the bases of the Feynman formalisms, assuming electromagnetic interaction acting 

•rH ■ at each vertex of the ladder diagram Fig.(l). Expanding the transition matrix T of the (ep) scattering in terms of 

y\ ' transition ladder diagrams r(") so that. 



t=Y[t, (1) 



{c"} and {r(")} are the coefRcients of expansion and transition matrices of the adder diagram of order n. T^") has 
the form, 

Where the factor stands for the Green's propagator of the virtual intermediate state number j. Vjj-i is 

the transition probability from the state j — 1 to j, so that, Vjj^i =< [F > . The screening Coulomb field 
working at each vertex has the form, V = —Uq ^ so that, each internal integration, Vkk' has the form; 

< k\V\k' ~Uo[2TT^{a^ + \k~ k'\^)]. (3) 
Hence, a one step ladder diagram has a transition matrix T'^^\ 



1 



TW = -C/o[27r2(a2 + |fci-fc^|2)] (4) 

ki and kf are the momentum of the initial and final states. Similarly, the two step and the three step ladders will 
have the forms, 

rp{2) _ _'yTr2Tj2 f (r\ 

°J (fci - fc2 - ie)(a2 + \ki - k[\^){a^ + \ki - k'J^) ^ ' 



y(3) 
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fc2 - ie)(a2 + \ki - k[\^){a^ + \ki - /e}|2) 

(^oV47r^)d fci d k2 

{k-kl- ie){a'^ + \k2 - fc;|2)(a2 + \k2 - k'f\'^){k - kj - ie)' 

(/e-fc?-k)(a2 + |fc,-fc;|2) 

The integrals in Eqs.(5;6) are carried out by the Dalitz integrals [12]. The expansion coefficients {c,} in Eq.(l) 
are determined by fitting with the experimental data to estimate the statistical weight factor of each ladder diagram 
contributing the reaction. According to the usual formalism for the dynamics of the particle reactions [13] , the phase 
space integral is an integration of the square of the transition matrix 

T{{ki}) over a set of allowed values {h}. 

^"(«)= J E^S'ik, + kb-Y.''i)\T{h}\'' (7) 

i i—1 

ka, kb are the momenta of the interacting particles in the initial state, while {fcijbeing the momenta of the particles 
in the final state. For the problem under consideration of the elastic scattering, where only two particles in the final 
state n = 2, with the normalization convention, the reaction cross section for ka + ki, ^ ki + k2 is 

A is a standard function defined as,A(.T, y, z) = {x — y — — Ayz . One may immediately write the differential 
cross section in the center of mass system CMS 



C^O- 1 k-Y 2 



d^ 647r2s ka 
It is more convenient to use the invariant cross section, 

da 1 ,„,, 



T|2 (9) 



Q2 (10) 



dt 167rA(s, m2, m2) ' 
Eqs. (9), (10) as well as Eq.(l) are in a relevant form to compare with the experimental data. 



III. THE MULTI PERIPHERAL MODEL (MPM) 

In this section, we shall deal with the problem of the particle production in the (ep) inelastic scattering in view of 
the multi peripheral collision [14]. A factorizable transition matrix T is assumed in the form, 

ra-l 

r=nri (11) 

i 

Each particle in the final state is produced at a specific peripheral surface as shown by the Feynman diagram Fig. (2) 
with electromagnetic transition matrix Tj, may be written in a suitable parametric form, 

r=^— (12) 

Oli + ti 
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where tj is the four vector momentum transfer square at the i peripheral surface, ai is the electromagnetic 
peripheral parameter characterizing the surface number i, and determine to conserve the total energy. The advantage 
of this technique is to reduce the many body problem into (n — 1) iterative diagrams, each of them has only two 
particles in the final state. For example, the i*'* diagram has two particles in the final state, the first one is the particle 
number {i + 1), and the other one has an effective mass Mj, equivalent to the rest of the i-particles of the system. 
The square of the 4- vector momentum transfer ti is kinematically calculated as, 

U = {P^^ -Pi - ■■■■-Pi) 



ti = m^+i +Mf- 2i;Wfc? - 2P^^Ki cos 0, (13) 

ruj+i is the rest mass of the particle number {i + 1) produced in the i*'' iteration. Ki and are the 3-vector 
momentum and the total energy of the effective mass Mi. and are the corresponding figures for the leading particle 

acting at the i*'' peripheral surface. The recursion relation of Pi^'is given by, 

Pi'^=X'/\Mf,ti,ml)/2Mi (14) 

The leading particle for the first peripheral surface is given by 

P„ = Pi") = AV2(s, ml, ml)/2^s (15) 

The multi peripheral parameters play important role in converging the particles in phase space and conse- 
quently, control the energy of the particles in final state. So that the values of are adjusted to conserve the total 
energy. The energy Ei of the particle number i is related to its rapidity t/j through the relation, 

Ei = nit cosh(yj) 



mt = ^JP^+mi (16) 
so that the total energy of the particles in the final state is 

C = — [mt cosh{y){da/dy)dy (17) 

which is a function of the parameters {ai} should be compared with the total center of mass energy ^/i of the 
initial state. We first start with n = 2 to get ai, which is inserted again in the case n = 3 to get a2 and so on. 

These are repeated sequentially to get the values of the rest parameters up to oin-i- The values {cti} depend on the 
particle multiplicity n rather than the energy y/s. The phase space integral /„(s) is then calculated as in Eq.(7) after 
transforming the integral variables from ki to ti, 



2M„ f\AP).''Ju, Jt7 . ti-i+ai-i^ 



H^ = J2mj (19) 

tf are the upper and lower limits of ti corresponding to cos{Oi) = ±1. 

tf=ml + ml - 2EaEb ± 2P,Pt (20) 

Multiple integrals in Eq.(18) are carried out using a Monte Carlo program, through which all possible distributions 
of the physical quantities are easily found. 
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IV. MONTE CARLO 



A Monte Carlo program GENE2 [15] is developed to simulate events according to the multi peripheral diagrams. It 
includes 3-gcncrators. The generator Gn{s) for the multiplicity of particles in the final state, the generator GM{s,n) 
for the invariant masses produced at the n—1 possible peripheral surfaces and finally the dynamic generator Gt{s, n, M) 
which generates the values of the momentum transfer square ti of the i*^ surface according to electromagnetic transition 
matrix. The program GENE2 executes the 3- generators in ordered sequences, and simulates all the kinematical 
variables of the n- particles in the final state. 

A general form of the algorithm used to generate an event x with probability f{x) is, 

r = / f{x)dx/ [ f{x)dx (21) 

Where a and h arc the boundaries of the physical region at which x is defined, r is the default random number with 
uniform distribution between the limits {0, 1}. 

A. The multiplicity generator Gn(s) 
It is assumed that the multiplicity distribution P{n) has a Gaussian form [16] , 

P(n) = ^exp[-(n-n)Va^] (22) 



where n is the average multiplicity which depends on the center of mass energy and cr is the dispersion. 

n= a log(s) + b ,a = c n +d (23) 
the integration limits extend from zero to infinity. This makes Eq.(21) read; 

r = Erf (n- n) (24) 
Erf{x) is the error function. The solution n of Eq. (24) defines the multiplicity generator of the reaction. 

B. The invariant mass generator GM(s,n) 

According to the multi peripheral diagrams Fig. (2), it is assumed that the peripheral surface number i, produces 
the particle number i + 1 and the invariant mass Mj, equivalent to the eff'ective mass of the system of the rest of the 
i— particles. The value of Mi satisfies the relation, 

J2^i<Mi<Mi+i (25) 

so that the algorithm of the generator GM(s,n) is, 

i i 

Mi = ^ ruj + r(Mi+i -^nij) (26) 

3 J = l 

For a system of n-particle final state, we start with M„ = ^/s, then generate the value of M„_i. Eq.(26) is the 
recursion relation to generate all values of Mj, i = n — 1, 2. 



C. The dynamic generator Gt(s,n,M) 

The algorithm used to generate the values of the square of the momentum transfer ti should simulate a probability 
distribution which is proportional to the square of the transition matrix Tj defined by Eq.(12). Inserting this in 
Eq.(21), we get the U generator as, 

ti = {r[{tt + - {t- + ai)-^ + (i+ + ai)-^}-^ - at (27) 
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V. RESULTS AND DISCUSSION 



A. The elastic scattering 

The problem of elastic scattering is treated in this article using the ladder diagrams of multi steps. The transition 
matrix T^"^ is calculated for ladder diagrams forn = 1, 2 and 3 using the Dalitz and Feynman integrals. The results are 
demonstrated in Fig. (3). The potential field acting at each vertex of the diagram is assumed to be of electromagnetic 
nature with screening factor due to the pion current inside the proton target. The screening effect limits the infinite 
range of the coulomb potential. It is found that the terms of the transition matrix T^"^ form a converging series. The 
first term of which decreases rapidly with the momentum transfer square t = —Q^. while the higher order terms are 
slowly varying functions. The asymptotic behavior of T*^"-* follows a power law at extremely high energy. 

The expansion coefficients {c;} are determined by the fitting method with the SLAC experimental data [4] in the 
range t ^ 1 — 3{GeV/c)'^. The results are shown in Fig. (4). 

It is found that only two terms of the series are sufficient to represent the reaction. The maximum probable 
diagram corresponds to n = 1. On the other hand, the SLAC-E-136 experiment, Fig. (5) in the momentum transfer 
range t 3 — 31{GeV/c)^ needs more terms of the series, with a most probable diagram corresponds also to n = 1. 
The values of the expansion coefficients are given in Table (1). 

Table(l) The branching ratios of multi photon exchanged in ep collisions. 



Reaction 


1-Step 


2_Steps 


3-Steps 


Rcf.[3] 


91 


9 




SLAC-E-136 


86 


12 
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Figure (4) and (5) show good agreement between the experimental data and the prediction of the multi photon 
exchange mechanism. It is clear that the reactions of high momentum transfer needs more terms of ladder diagrams. 

B. The inelastic scattering. 

Here we use the iterative multi peripheral diagram as in Fig.(2-a,b). The advantage of this is to simulate a complete 
inclusive reaction. We used the Monte Carlo program GENE2 to generate the n-particles in the final state. The multi 
peripheral parameters ai, i = 1,2, 1 are determined to conserve the total center of mass energy y/s. The values 
of the parameters {ai} depend not only on ^/s but also on the degree of pcripherality i. The first few peripheral 
surfaces possess parameter values which confine the particle production in very narrow cone. The production cone 
angle gets wider for higher order surfaces. The values of are displayed in Fig. (6) as a function of the multiplicity 
n, of the number of particle in the final state for electron lab energies 50, 100 and 2GGG'ey. In all cases the value of 
decreases slowly with n, up to a critical multiplicity nc after which a sudden drop is observed. A parametric relation 
is obtained for ai as a function of the electron lab energy E and n as, 

a„_i = (0.226 - 0M5E + 0.0003^^)^ ^ o.441 log(^) - 2.29, n < ric. 
a„_i = {0mi9E - 0.8685) exp[(-.000688£; -|- 0.327276)n], n > ric 

Uc = 0.4846£; + 1.0769 (28) 

The critical value ric corresponds to the peripheral surface at which enough energy is transferred, that is sufficient 
to make phase transition from the nuclear matter to the quark gluon state. The model is applied to the data of the 
experiment SLAC-E-133 corresponding to electron lab energies 9.744,12.505,15.730,18.476 and 20.999Gey. Fig.(7) 
shows the rapidity distributions for particles produced at 9.744GeV electron lab energy as calculated by (MPM) at 
peripheral surfaces n = 2, 3, 5 and 10. The result for n = 2 (two particle final state) shows two non symmetric clear 
peaks corresponding to forward and backward emission. The electron peak is relatively narrower than the proton one 

. This case represents the most peripheral collision with low momentum transfer. As the multiplicity increases the 
two peaks get broader and interfere through each other. They completely interfere and form only one peak at the case 
of high multiplicity n = 10 corresponding to high momentum transfer. This represents the most central collision, i.e. 
collision between the electron and the proton core. The missing mass of the recoil nucleon is also calculated through 
the Monte Carlo program GENE2, according to the relation. 
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W'^ = M^ + 2Mp{E - E') - AEE' sin^ - 



(29) 



The differential cross section da/dCl as a function of is compared in Fig. (8) with the experimental data of 
energies 9.744, 12.505 and 15.730GeV. Fair agreement is obtained for the first reaction only. The deviation increases 
as the electron energy increases. This may be due to the fact that the model in hand has ignored the relative motion 
of the core inside the nuclcon target. This relative motion let the recoil nucleon not has a unique value of missing 
mass, but instead, it posses a Gaussian -like distribution around W"^ = M^. The peak is wider at higher t. 



VI. CONCLUSION 



1- Data of elastic scattering may be reproduced in a wide range of the energy transferred using the (MPEM) with 
an expandable transition matrix representing ladder diagrams. 

2- A transition matrix of order n, may be approximated as the power form of the one step ladder diagram at 

extremely high energy. 

3- The power series of the transition matrix contains number of terms increasing with the energy transferred of the 
reaction. 

4- Inelastic ep scattering are successfully described by the multi peripheral model. 

5- The peripheral parameters are the dynamic parameters, which control the convergence of the particles in phase 
space. A critical value of peripheral parameter dcterinines the critical surface at which phase transition may occur 
from the nuclear matter density to the quark gluon phase. 

6- The rapidity distribution of the particles in the final state shows two peaks representing the forward and backward 
production. The peaks interfere together as the multiplicity increases in the final state. 
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